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ABSTRACT 


Multispectral imaging (MSI) data collected at multiple angles over shallow water provide 
analysts with a unique perspective of bathymetry in coastal areas. Observations taken by 
DigitalGlobe’s WorldView-2 (WV-2) sensor acquired at 39 different view angles on 30 
July 2011 were used to determine the effect of acquisition angle on derived depth. The 
site used for this study was on the island of Oahu, focused on Kailua Bay (on the 
windward side of the island). Satellite azimuth and elevation for these data ranged from 
18.8 to 185.8 degrees and 24.9 (forward-looking) to 24.5 (backward-looking) degrees 
(respectively) with 90 degrees representing a nadir view. Bathymetry were derived 
directly from the WV-2 radiance data using a band ratio approach. Comparison of results 
to LiDAR-derived bathymetry showed that varying view angle impact the quality of the 
inferred bathymetry. Derived and reference bathymetry have a higher correlation as 
images are acquired closer to nadir. The band combination utilized for depth derivation 
also has an effect on derived bathymetry. Four band combinations were compared, and 
the Blue & Green combination provided the best results. 
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I. INTRODUCTION 

A. PURPOSE OF RESEARCH 

Imagery data acquired from satellites are widely useful in the field of 
oceanography. These data contribute to a better understanding of the dynamics of ocean 
circulation, are useful for monitoring climate change, can be used for navigation and 
fisheries management, and are also helpful for improving models of ocean circulation, 
air-sea interaction, weather forecasting, and climate [Sanford el al., 2011]. 

Information about shallow water bathymetry is beneficial to scientists or groups 
that require knowledge of ocean depths in a particular coastal location. The use of 
multispectral imagery (MSI) data has been shown to adequately determine depths of 
remote coastal areas, as SOund Navigation And Ranging (SONAR) or Light Detection 
And Ranging (LiDAR) soundings or other bathymetric data may not be available prior to 
arrival. 

There are many benefits to using remote sensing data acquired from satellites. 
Spaceborne sensors can collect spectral data over extremely large areas which can be 
advantageous to regions not accessible on foot, or for areas denied by hostile forces. 
These data can also be acquired much more rapidly due to the number of space platforms 
collecting on a regular basis. Finally, improvements in collection capabilities steadily 
grow as space platform production increases and sensor technologies evolve. This 
capability is extremely valuable to the Navy and Marine Corps, who routinely conduct 
and plan worldwide amphibious operations. Coastlines are highly dynamic environments 
and mission success is dependent upon knowing the underwater terrain. Being able to 
acquire bathymetry in a contested or denied region greatly improves the probability of 
success for these amphibious operations. 

The purpose of this research was to integrate the use of multiple satellite image 
acquisition angles over one location, and detennine what role these varying angles play in 
bathymetric depth determination. Analyses of 39 WorldView-2 (WV-2) images acquired 
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over Kailua Bay on the windward side of Oahu, Hawaii were used to reach a conclusion 
about the accuracy of bathymetric derivation from MSI. 

B. SPECIFIC OBJECTIVES 

The objective of this study was to test the potential of bathymetric derivation 
using WV-2 imagery acquired at multiple angles, and then report upon the role that 
image acquisition angle plays in depth detennination. 

The motivation for this work originated from the need to determine bathymetry 
from only a single spectral image of the coastal region in question. There is no guarantee 
that this image will have been acquired at optimal viewing geometry, e.g., nadir. Without 
sufficient time to task a satellite and acquire data at a particular viewing angle before a 
site visit, results from this research strive to provide a better understanding of how to 
manipulate image data in order to obtain a better understanding of bathymetry. Analysis 
of multi-angle MSI data was used to quantify the effects of varying satellite acquisition 
angle on accurate determination of bathymetry. 
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II. BACKGROUND 


This research focused on the use of passive, optical remote sensing systems for 
image acquisition. Passive sensors are those that rely on incoming solar radiation to 
illuminate the targets on land and in the water [ Camacho , 2006]. Optical sensors are 
those that focus on the visible (VIS), near infrared (NIR), and shortwave infrared (SWIR) 
portions of the electromagnetic (EM) spectrum to observe radiation from targets 
[ Camacho , 2006]. The following sections will explain these principles in more detail, as 
well as further discuss some issues that affect remote sensing of ocean environments, as 
well as imagery acquired at multiple angles. 

A. PRINCIPLES OF RADIATIVE TRANSFER 

A remote sensing instrument or sensor receives energy that is reflected from the 
surface of the Earth. This energy is affected by interactions of light with the atmosphere 
and water, as well as any particulate matter in the water column [ Camacho , 2006]. It is 
imperative to understand the basics about the EM spectrum, radiation, the variety of 
possible interactions, and line spectra. 

1. Electromagnetic Spectrum 

The EM spectrum is a conglomeration of a number of classifiable spectral 
regions. These include: gamma rays, X-rays, ultraviolet (UV) light, VIS light, IR light, 
microwaves, and radio waves (Figure 1) [Olsen, 2007]. 
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Visible Portion of the electromagnetic Spectrum from 400 to 700 nm 



Figure 1. The Electromagnetic Spectrum (From Jensen [2007]) 

This research focuses on the VIS and NIR portions of the spectrum, extending 
from approximately 400 nm to 750 nm, and 750 nm to 1000 nm (respectively). 

2. Spectral Signatures 

Single atoms or molecules emit light in the form of line spectra, and an atom that 
is well isolated will radiate a discrete set of frequencies called a line spectrum. The 
wavelengths radiated and absorbed are specific to that atom or molecule, and are 
presented as spectral signatures. These signatures are used to determine the composition 
of radiating or absorbing gases and other materials [Olsen, 2007]. 
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Spectral features, such as those shown in Figures 2 and 3, enable analysts to better 
differentiate between materials on the ground. Based on specific values, peaks, and 
troughs within the graph, one can determine the basic makeup of a material within the 
image. 


A (%) 



Figure 2. Hyperspectral Imagery (HSI) provides more complete spectral signatures than 
MSI because there are a much larger number of bands. Albedo values for 
different forms of algae and coral are shown (From Maritorena et al. [1994]; 

Camacho [2006]) 
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Figure 3. MSI has spectral signatures that look less complete compared to HSI because 
there are fewer bands (spectra are taken from a scene used in this research and are 

in nanometers) 


3. The Four Fundamental Energy Interactions with Matter 

Electromagnetic radiation may be transmitted, reflected, scattered, or absorbed. 
The proportions to which these interactions occur depend on the compositional and 
physical properties of the medium, the wavelength or frequency of the incident radiation, 
and the angle at which the incident radiation strikes a surface [Avery and Berlin, 1992; 
Olsen, 2007]. 


a. Transmission 

With transmission, shown in Figure 4, incident radiation passes through 
matter without measurable attenuation. Different densities of the material, however, can 
cause radiation to be refracted or deflected from a straight-line path, and will also alter 
the velocity and wavelength. The change in EM radiation velocity can be further 
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explained by the index of refraction, which is the ratio between the velocity of the EM 
radiation in a vacuum (a perfectly transparent medium) and its velocity in a material 
medium. 


The index of refraction for a vacuum is equal to 1, and can never be less 
than 1 for any substance [Avery and Berlin, 1992]. 



Figure 4. Transmission (After Avery and Berlin [1992]) 
b. Reflection 

Reflection, or specular reflection, occurs when surfaces are smooth 
relative to the wavelengths of incident radiation (Figure 5). Specular reflection, the 
process by which incident radiation bounces off the surface of a substance in a single and 
predictable direction, causes no change to either the EM radiation velocity or wavelength 
[Avery and Berlin, 1992]. 
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Figure 5. Reflection (After Avery and Berlin [1992]) 
c. Scattering 

Scattering, or diffuse reflection, takes place when incident radiation is 
dispersed in unpredictable directions. Scattering occurs when surfaces are rougher 
relative to the wavelengths of incident radiation, as shown in Figure 6. The velocity and 
wavelength of EM waves, however, are not affected by scattering [Avery and Berlin, 
1992]. 



Figure 6. Scattering (After Avery and Berlin [1992]) 
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(1) Bidirectional Reflectance Distribution Function (BRDF). The 
Bidirectional Reflectance Distribution Function (BRDF) describes the scattering 
characteristics of a material by describing how light from a source, incident on a target, is 
reflected in a given direction. In order to calculate BRDF, one looks at the ratio of the 
energy scattered by the target in a particular direction, dependent on wavelength, over the 
energy that was incident on the target from a particular direction. All combinations of 
inbound and outbound energy directions are then integrated. BRDF is dependent on both 
the angle of incidence and on the angle of reflectance. Angles are defined in a Cartesian 
coordinate system by a polar angle, 9, measured from the surface normal, as well as an 
azimuthal angle, (p, measured from the x-axis (Figure 7) [McConnon, 2010], 



y 


Figure 7. Reflection geometry used in the definition of BRDF 

(From McConnon [2010]) 


d. Absorption 

Absorption is the process by which incident radiation is taken in by a 
medium. This occurs when a substance is opaque to the incident radiation, a portion of 
which is converted to internal heat energy, then emitted or reradiated at longer thermal 
infrared wavelengths (Figure 8) [Avery and Berlin, 1992]. 
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Emission 



Figure 8. Absorption (After Avery and Berlin [1992]) 

B. INTERACTIONS OF LIGHT WITH THE ATMOSPHERE 

There are three limiting factors that the atmosphere introduces into the field of 
remote sensing. These include atmospheric absorption, scattering, and turbulence [Olsen, 
2007]. Atmospheric turbulence, however, will not be discussed here because its impact is 
greater for telescopes looking up from the earth through the atmosphere than for sensors 
looking down. Also, it is important to keep in mind that, overall, the atmosphere is more 
transparent in the long-wave IR (1100 to 1200 nm) than in the VIS spectrum (400 to 750 
nm) [Olsen, 2007]. 

1. Atmospheric Absorption 

Atmospheric absorption is mainly dependent on wavelength, and is most affected 
by water (between 500 and 700 nm), carbon dioxide, and ozone (near 1000 nm) [Olsen, 
2007]. Even after sensor calibration has occurred, the shape of the measured spectra 
relative to the underlying reflectance spectra are altered because of atmospheric 
absorption lines [Eismann, 2012]. Common absorption spectra can be seen in Figure 9. 
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Figure 9. The first four graphs show the absorption characteristics of N 2 O, 

0 2 and 0 3 , C0 2 , and H 2 0, while the bottom graphic depicts the cumulative result 
of all these constituents being in the atmosphere at one time (From Jensen [2007]) 

2. Atmospheric Scattering 

Atmospheric scattering is primarily caused by collisions between aerosols and 
particulates like dust, fog, and smoke in the atmosphere [Olsen, 2007]. There are two 
main consequences of atmospheric scattering: reduction of radiant energy and unwanted 
gain in the sensor [ Martin , 2004; Camacho, 2006]. Scattering can also be divided into 
two types: Rayleigh and Mie scattering. Rayleigh, or molecular, scattering is primarily 
caused by oxygen and nitrogen molecules, i.e., particles whose effective diameters are 
less than the wavelengths of interest. As the size of the particle increases, the scattering 
processes move towards Mie scattering. This type of scattering occurs when the diameter 
of particles is similar to the wavelengths of the energy being radiated [ Camacho, 2006; 
Olsen, 2007]. 


11 



3. 


Atmospheric Compensation for Spectral Imagery 


The effects the atmosphere has on incoming radiation reaching a remote sensing 
sensor results from a variety of processes. Five possibilities are shown in Figure 10 [Kay 
et al., 2009]. 



— ► A Single or multiple scattering in the atmosphere, by molecules or aerosols 

— ► it Scattering from the atmosphere to the water surface and then reflection to the detector often 
termed “sky glint" 

—► C Reflection from whitecaps on the sea surface 

—► T) Specular reflection from water surface, with direct transmission through the atmosphere from 
the sun to the surface and from the surface to the detector - this is termed sun glint' in the 
context of this review. 

— ► F. Transmission through the atmosphere and air-water interface followed by scattering or 
reflection below the water surface and transmission back through the atmosphere to the 
detector. 


Figure 10. Diagram showing routes by which light can reach a remote 

sensing detector (From Kay et al. [2009]) 


Kay et al. [2009] assumes that if these five processes are predominantly 
responsible for the sensor-received signal, then: 


( 1 ) 



'water 


where T is the transmittance of the atmosphere along the sensor view direction. The path 
L atm is the radiance that arrives at the sensor via atmospheric paths (single or multiple 
scattering in the atmosphere, by molecules or aerosols). Other terms L s k y , Twhitecap, Tgiint, 
and Lwater are the radiances just above the water surface for light traveling by sky glint, 
whitecap, sun glint, and water-leaving routes. These terms depend on the wavelength, 
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along with other factors. The term L walcr contains the infonnation about water column 
and benthic features, and needs to be separated from the other terms if this infonnation is 
going to be retrieved [Kay et al., 2009]. 

C. INTERACTIONS OF LIGHT AND WATER 

Optical oceanography is vital for addressing problems such as photosynthesis, 
ecosystem dynamics, ocean health, seawater clarity, underwater imaging, biogeochemical 
cycling, carbon budgets, upper-ocean thennodynamics, and climate change [Dickey et 
al, 2011], 

The simplest optical interactions occur at the boundary of the atmosphere and 
ocean, and are governed by Snell’s law and the Fresnel equations. In order to accurately 
model the very complex interactions between water molecules and other constituents, 
radiative-transfer methods and precise measurements are required [Dickey et al., 2011]. 

Optical remote sensing is an important tool utilized for monitoring marine 
environments. Since light is readily absorbed by water, remote sensing is most 
successful in places with very shallow, clear water, usually up to depths of 30 or possibly 
40 meters. Sea water contains dissolved and particulate matter, with varying 
concentrations (both spatially and temporally) throughout the water column [Mobley, 
1994; Camacho, 2006]. Optical properties of the water column can be divided into two 
classes—inherent optical properties (IOPs) and apparent optical properties (AOPs) [Smith 
and Baker, 1981; Camacho, 2006]. These classes are further discussed in sections to 
follow. 


1. Inherent Optical Properties (IOPs) 

An optical property is inherent if it depends only upon the medium and is 
independent of the ambient light field within that medium [Mobley, 1994; Camacho, 
2006]. As light enters the water column, it interacts with particles and sediment, which 
cause the incident light to be altered by scattering or absorption [Thomas and Stamnes, 
1999; Camacho, 2006]. It is these scattering and absorption characteristics that define 
IOPs of water [Camacho, 2006]. 
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2. Apparent Optical Properties (AOPs) 

An optical property is apparent if it is dependent on the medium, as well as the 
directional structure of the ambient light field. Similar to IOPs, AOPs are also dependent 
on the dissolved particles and sediment in the water column. Unlike IOPs, however, 
these properties cannot be measured in situ because they depend on the ambient radiance 
[Mobley, 1994; Camacho, 2006], 

3. Sun Glint Correction for Spectral Imagery 

Sun glint is the specular reflection of light directly transmitted from the upper side 
of the air-water interface [Kay et al., 2009]. It typically fonns bands of white along wave 
edges on the windward side of nearshore environments [Heclley et al., 2005], As it is a 
serious confounding factor for remote sensing of water column properties and benthos, 
researchers have spent a good deal of time creating and testing techniques to estimate and 
remove the glint radiance component from imagery. Sun glint, which is a function of sea 
surface state, sun position, and viewing angle, occurs in imagery when the water surface 
orientation is such that the sun is directly reflected toward the sensor. The component of 
sensor-received radiance can be so high that it becomes impossible to retrieve any 
information about the ocean environment [Kay et al., 2009]. Examples of sun glint 
appearing in imagery are shown in Figures 11 and 12. 
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NOAA GOES 12 


091202 0845 UTC 


NASA GSFC GOES Project 



Figure 11. An example of sun glint on the ocean’s surface (image captured by the 

Geostationary Operational Environmental Satellite (GOES) on 2 December 2009) 
(Image from http://rsd.gsfc.nasa.gov/goes/text/hotstuff.html) 



Summer Sunglint in the Sargasso Sea 
dawns on Cuba and Florida 
1200 GMT on 22 June 2000 
GOES-West visible image 


GOES project NASA-GSFC 


i% Q 


Figure 12. Another example demonstrating the effect of sun glint on the imager 
(image captured by the GOES satellite on 22 June 2000) (Image from 
http ://rsd. gsfc .nasa. gov/goes/text/ho tstuff.html) 
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Sun glint causes problems within the imagery domain. A case study performed 
by Goodman et al. [2008] reported that uncorrected glint in high resolution imagery led 
to errors as large as 30% when measuring ocean depth. For airborne surveys, optimal 
flight paths and directions can be chosen based on the time of day, but this is more 
difficult in the case of satellite imagery [ Goodman et al., 2008]. 

Glint correction methods have been developed to improve image accuracy for two 
main categories of water types: open ocean imaging, and higher resolution coastal and 
aerial applications. The aim for both cases is the same—to estimate the glint contribution 
to the radiance reaching the sensor, and then subtract it from the received signal [Kay et 
al., 2009], 

The set of methods utilized for coastal images is meant to be used with pixel sizes 
of less than approximately 10 meters. Data from the NIR band are used as an indication 
of the amount of glint in the received signal. The spectrum from a section of the image 
containing deep water is then used to establish the relationship between the NIR and glint 
radiances [Hochberg et al., 2003; Hedley et al., 2005; Lyzenga et al., 2006; Goodman et 
al., 2008; Kay et al., 2009]. The assumption of a linear relationship between NIR 
brightness and the amount of sun glint in the VIS bands holds because the real index of 
refraction is approximately equal for NIR and VIS wavelengths [Mobley, 1994; Hedley et 
al., 2005]. 

When no sun glint is present, the radiance received by a satellite-borne sensor is 
dominated by atmospheric scattering. Light from paths through the air constitute over 
80% of received radiance, water-leaving paths make up approximately 15%, and 
reflected light only 1 to 2% [Sturm, 1981; Kay et al., 2009]. Sun glint can increase the 
reflected radiance by a factor of 2 or more, and the worst case glint can saturate the 
sensor, making it impossible to retrieve the water-leaving radiance for those pixels [Kay 
et al., 2009]. 
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D. MULTI-ANGLE RESEARCH 


Multi-angle remote sensing capabilities offer a number of advantages with respect 
to a single shot dataset. Multi-angular data fusion has been shown to allow: 

• The exploitation/investigation of BRDF, 

• The extraction of digital height maps (DHMs), 

• Atmospheric parameter retrieval, 

• Classification improvement, etc. [Pacifici et al., 2011] 

The implications of adding height data and multi-angle MS reflectance, both 
derived from the multi-angle sequence, to the textural, morphological, and spectral 
information of a single MS image has been investigated by Longbotham et al. [2011], It 
was detennined that multi-angle collections significantly increase the dimensionality of 
the data available for a single target, as well as allow for differentiation of classes not 
typically well identified by a single image [Longbotham, Chaapel et al., 2011]. 

For the particular Longbotham et al. [2011] study referred to here, the analyzed 
data sequence was collected over urban areas of Atlanta, GA in December 2009. The full 
dataset contained 27 images, with satellite elevation (the angle between the horizon and 
satellite as viewed from the image target) as low as 25 degrees, but the study only utilized 
those images with a moderate (relative to the satellite’s capability) off-nadir observation 
angle of less than 30 degrees. This left the group with 13 images between 57 degrees 
(forward-looking) to 81.5 degrees (most nadir) to 59 degrees (backward-looking). 
Results showed a 27% improvement in classification accuracy for spatial experiments, 
and also a 14% improvement in classification accuracy for spectral experiments. Figure 
13 shows the ground observed azimuth and elevation of the WV-2 satellite for each 
image observation, as well as the azimuth and elevation of the sun during the multi-angle 
sequence used for the research [Longbotham, Chaapel et al., 2011]. Further sections will 
discuss why WV-2 is capable of acquiring multi-angle imagery. 
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0° (N) 



180° 


Figure 13. Ground observed azimuth and elevation of the WV-2 satellite 

(black markers) for each image observation as well as the azimuth and elevation 
of the sun (yellow circle) during the 2009 Atlanta, GA multi-angle sequence 
acquisition. Azimuth is plotted angularly clockwise (North: 0 degrees, East: 90 
degrees, South: 180 degrees, West: 270 degrees), and elevation is plotted radially 
from the center (ground nadir: 90 degrees, ground horizon: 0 degrees) (From 
Longbotham, Chaapel et al. [2011]) 


There is significant improvement shown over the baseline classification when 
using a multi-angle WV-2 sequence. In terms of spatial classification, Longbotham et al. 
[2011] successfully demonstrated the ability to differentiate between classes like bridges 
and man-made structures, which are generally difficult to classify because they are 
spectrally similar to ground-level classes of the same material [ Longbotham et al., 2011]. 
Improvements were also made in terms of spectral classification—the group was able to 
distinguish between classes that may prove to be valuable in land-use classification, such 
as moving versus parked vehicles. The ability to differentiate between spectrally similar 
classes such as paved parking lots and paved highways were also made possible using the 
multi-angle sequence [ Longbotham, Bleiler et al., 2011]. 
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E. PREVIOUS WORK AT THE NAVAL POSTGRADUATE SCHOOL 
1. Work in Satellite Derived Bathymetry 

a. Depth Analysis of Midway Atoll Using QuickBird Multi-Spectral 
Imaging over Variable Substrates (Camacho, 2006) 

Camacho [2006] used QuickBird MSI over Midway Atoll, Northwestern 
Hawaiian Islands in an attempt to analyze depth and identify variable bottom-types in 
shallow water. The objective of this study was to use an MS image to categorize benthic 
substrates based on spectral characteristics and ground truth data collected in situ. With 
this image, a ratio of reflectances (the “ratio method”) was used to extract depth 
separately over variable substrates. The motivation for this work originated from a 
limitation pointed out by another Naval Postgraduate School thesis student, Clark [2005], 
whose results demonstrated that the ratio method proved to be sensitive to bottom type. 
The method produced shallower depths over bottom types with low albedo and deeper 
depths over bottom types with high albedo. Clark [2005] also noted that sun glint had an 
effect on the overall results. The maximum depth that could be obtained using the ratio 
method was 15 meters [Clark, 2005]. 

Camacho [2006] utilized methods by Lyzenga [1978] and Stumpf et al. 
[2003] to extract water depth and bottom type information, and also used the ratio 
method to retrieve accurate depths over variable bottom types and low-albedo 
environments (respectively) [Lyzenga, 1978; Stumpf et al., 2003; Camacho, 2006], 

b. Depth Derivation from the WorldView-2 Satellite Using 
Hyperspectral Imagery (Loomis, 2009) 

Loomis [2009] created a simulation before WV-2’s launch to determine 
the usefulness of the Yellow band for depth derivation. Loomis [2009] utilized data from 
the Airborne Visible/InfraRed Imaging Spectrometer (AVIRIS) sensor over Kaneohe 
Bay, Hawaii, and then processed the imagery using the Stumpf et al. [2003] ratio method 
to detennine bathymetry. The Green & Blue, Yellow & Green, and Yellow & Blue band 
combinations were compared to ground truth bathymetry derived from a digital nautical 
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chart. Results indicated that use of the Yellow band improved the accuracy of derived 
depths, especially in shallow water [Loomis, 2009]. 

c. Contributions to Remote Sensing of Shallow Water Depth with 
the WorldView-2 Yellow Band (Madden, 2011) 

Following Loomis’ [2009] work, Madden [2011] analyzed bathymetry in 
Tampa Bay, Florida’s shallow water using WV-2 imagery. The Yellow band was 
combined separately with the Blue, Green, and Red bands, and then compared to the 
more traditional Blue & Green and Green & Red combinations. Madden’s [2011] results 
showed that the addition of the Yellow band provided more information about 
bathymetry, but less sensitivity to bottom type in two of three transect lines used 
[Madden, 2011]. 

2. Work in Multi-Observation Imagery 

a. Coastal Bathymetry Using 8-Color Multispectral Satellite 
Observation of Wave Motion (McCarthy, 2010) 

McCarthy [2010] measured coastal bathymetry near Camp Pendleton in 
California by using wave motion as observed by WV-2. After principal component 
transforms were performed, successive change detection images contained both spatial 
and temporal information. Wave celerity could be detennined and depth inversion was 
performed. McCarthy [2010] measured the wavelength of a wave above a known depth, 
and then used the wave period method to detennine depth for other waves in the 
propagation direction of that particular wave. McCarthy [2010] detennined that the 
spatial resolution for this method was higher and possibly more accurate than the 
reference bathymetry used, especially in the surf zone [McCarthy, 2010], 

b. Coastal Bathymetry Using Satellite Observation in Support of 
Intelligence Preparations of the Environment (Myrick, 2011) 

My rick [2011] followed McCarthy’s [2010] work, and calculated water 
depths using linear dispersion relationships for surface gravity waves. Myrick [2011] was 
able to establish depth profiles out to 1 kilometer offshore, as well as derive depths up to 
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15 meters. Comparisons with United States Geological Survey (USGS) bathymetric 
acquisitions from 2009 show agreement within 5% in the surf zone and 1% outside of the 
surf zone [ Myrick , 2011]. 

c. Automating Nearshore Bathymetry Extraction from Wave 
Motion in Satellite Optical Imagery (Mancini, 2012) 

Mancini [2012] extracted nearshore depths for Waimanalo Beach, Hawaii 
from WV-2 optical imagery by means of automated wave kinematics bathymetry (WKB). 
Two sets of three sequential images, acquired at approximately 10 second intervals, were 
used. Depths from close to shore out to about 20 meters deep were generated. 
Comparisons to the Scanning Hydrographic Operational Airborne Lidar Survey 
(SHOALS) LiDAR bathymetry values showed WKB depths were accurate to about 0.5 
meters, with R values of 90%, and were frequently in the range of 10% to 20% relative 
error for depths ranging from 2 to 16 meters [ Mancini , 2012], 

d. High Spatial Resolution Bidirectional Reflectance Retrieval 
Using Satellite Data (McConnon, 2010) 

McConnon [2010] analyzed 15 WV-2 images acquired over Duck, North 
Carolina and 10 images of Pendleton, California to extract BRDF. Separation angles 
between the inbound, solar zenith angle, and the outbound reflectance angle were 
calculated. BRDF shifts across wavelengths and regions of interest (ROIs) were plotted 
against the previously mentioned values and then examined [McConnon, 2010]. 
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III. DATA AND METHODS 


A. DATA 

1. Satellite Sensor, Study Site, and Imagery Dataset 
a. WorldView-2 (WV-2) 

WorldView-2, launched on 8 October 2009 from Vandenberg Air Force 
Base, is DigitalGlobe’s third operational satellite after QuickBird and WorldView-1. The 
sensor is in a nearly circular, sun-synchronous orbit, and flies at an altitude of 
approximately 770 kilometers (Table 1). At nadir, the best possible spatial resolution is 
0.46 meters (panchromatic) and 1.84 meters (MS). What differentiates WV-2 from 
DigitalGlobe’s previous sensors is its 9 spectral bands—one panchromatic (PAN), 
ranging from approximately 450 to 800 nm (centered at 632 mn), and 8 MS, ranging 
from approximately 400 to 1050 nm (Figure 14, Table 2). These MS bands include: 
Coastal (centered at 427 nm), Blue (centered at 478 nm), Green (centered at 546 nm), 
Yellow (centered at 608 nm), Red (centered at 659 nm), Red Edge (centered at 724 nm), 
NIR-1 (centered at 831 nm), and NIR-2 (centered at 908 nm) (available at 
http://www.digitalglobe.com/about-us/content-collection#satellites&worldview-2). 
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Table 1. WV-2 Design and Specifications (After infonnation available at 
http://www.digitalglobe.com/about-us/content- 



collection#satellites&worldview-2) 

Launch Information 

Date: October 8, 2009 

Launch Site: Vandenberg AFB, CA, USA 

Orbit 

Altitude: 770 km 

Type: Sun synchronous 

Period: 100 minutes 

Sensor Resolution 

Panchromatic: 0.46 m GSD at nadir, 0.52 m GSD at 20 

degrees off-nadir 

Multispectral: 1.85 m GSD at nadir, 2.07 m GSD at 20 

degrees off-nadir 

Swath Width 

16.4 km at nadir 

Attitude Determination 

3-axis stabilized 

and Control 

Actuators: Control Moment Gyros (CMGs) 

Sensors: Star trackers, solid state IRU, GPS 

Pointing Accuracy and 

Accuracy: <500 m at image start and stop 

Knowledge 

Knowledge: Supports geolocation accuracy below 

Retargeting Agility 

Time to Slew 200 km: 10 seconds 

Revisit Frequency (at 

1.1 days at 1 m GSD or less 

40 degrees N Latitude) 

3.7 days at 20 degrees off-nadir or less (0.52 m GSD) 

Geolocation Accuracy 

Demonstrated <3.5 m CE90 without ground control 

Capacity 

1 million sq km per day 
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WV02 Relative Spectral Radiance Response 



Figure 14. WV-2’s relative spectral radiance response 
(From Updike and Comp [2010]) 

Table 2. WV-2 Sensor Bands (After information available at 
http://www.digitalglobe.com/about-us/content- 
collection#satellites&worldview-2) 


Band 

Band Width (nm) Center Wavelength (nm) 

Panchromatic 450 - 800 

632 

Coastal 

400-450 

427 

Blue 

450-510 

478 

Green 

510-580 

546 

Yellow 

585 -625 

608 

Red 

630 - 690 

659 

Red Edge 

705 - 745 

724 

NIR-1 

770 - 895 

831 

NIR-2 

860 - 1040 

908 
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The WV-2 system’s 8 MS bands are arranged in two arrays of 4 MS bands 
each (MSI and MS2). MSI includes Blue, Green, Red, and NIR-1. MS2 is comprised of 
Coastal, Yellow, Red Edge, and NIR-2. Imaging options are: PAN only, PAN + MSI, 
and PAN + 8 MS (MSI and MS2). There is a 0.3 second delay between MSI and MS2 
acquisitions, according to G. Miecznik (unpublished data, 2012). 

When combined, these bands are designed to improve the segmentation 
and classification of land and aquatic features beyond any other multispectral satellite 
imager. After WV-2’s launch, it was speculated that the increased agility and addition of 
the Coastal band would improve remote bathymetric measurements (mainly due to the 
Coastal band’s wavelength value, making it least absorptive by water). Analysts 
expected to be able to calculate depths up to 20 meters, and possibly even 30 meters 
using the Coastal, Blue, and Green bands. Once scientists were able to utilize WV-2 
imagery, the Coastal band has proven to be useful for the retrieval of water depth, true- 
color correction for human vision representation, chlorophyll absorption, and 
atmospheric scattering correction [Pacifici and Navular, 2011; Marchisio el al., 2011]. 
Also, with an average revisit time of 1.1 days and the ability to utilize 
off-nadir pointing, WV-2 was expected to better monitor ocean environments, which 
are highly dynamic and constantly changing (information available at 
http://www.satimagingcorp.com/media/pdf/Bathymetry_Datasheet.pdf). 

Bi-directional scanning is supported by WV-2. The onboard camera has a 
standard maximum look angle of approximately 40 degrees off-nadir and can slew across 
300 kilometers of the Earth’s surface in 9 seconds. This unique agility enables WV-2 to 
collect dense image sequences of a single target from several observation angles as it 
overflies an area [Longbotham el al., 2011]. 

b. Kailua Bay, Oahu, Hawaii 

This study focused on one location on the windward side of the Oahu 
coastline of Hawaii - Kailua Bay (Figure 15). The bay’s approximate latitude/longitude 
is: 21 degrees 24’ 29” N, 157 degrees 44’ 09” W. This particular beach was chosen 
because it was the least cloudy out of the 39 WV-2 images in the acquired dataset. 
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Figure 15. The Hawaiian Islands (left), focusing on Kailua Bay, Oahu (right) (From 

Google Earth) 


Kailua Bay is a carbonate reef-dominated embayment. There are two 
categories of benthic substrate found here: areas of carbonate sand and fossil reef 
hardgrounds, and reef habitats of coral and algae species. There is a sand-floored channel 
at the center of the bay, which cuts across the reef and connects the seaward and 
nearshore sand fields. Algae and corals grow on the plains [Isoun et al., 2003]. 


Areas with sand and fossil reef appear light-colored and are highly 
reflective in the WV-2 imagery. The coral and algae communities look dark and have 
low reflectance [Isoun et al., 2003]. Aerial images are shown in Figure 16. 



Figure 16. Oblique aerial imagery acquired over Kailua Bay on 10 November 2003 
on a Nikon Coolpix 5700 digital camera (settings and focus all automated) 
from a Cessna plane flying at about 2,000 feet (From 
http://www.soest.hawaii.edu/coasts/erosion/oahu/oblique.php) 
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c. 


Imagery Dataset 


A collection of 39 multispectral images of the windward side of Oahu, 
Hawaii was acquired by the WorldView-2 satellite on 30 July 2011 between 21:22:49Z 
and 21:28:54Z. Images were acquired at approximately 10 second intervals, covering 
just over 6 minutes (Figure 17). Section A of the Appendix contains a sample WV-2 
metadata file. 



Figure 17. Five examples out of the 39 WV-2 image acquisitions 

(labels are actual Image IDs) 


2. Bathymetry Data 

Bathymetry data of Kailua Bay were acquired from the University of Hawaii at 
Manoa, School of Ocean & Earth Science & Technology, Department of Geology and 
Geophysics (data are available at http://www.soest.hawaii.edu/coasts/ 
data/oahu/shoals.html and information about this dataset is available at 
http://www.soest.hawaii.edU/coasts/data/readme.html#shoals). According to the 
University of Hawaii website, these bathymetry data were collected from the SHOALS 
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website, as part of a survey conducted in 2000 (information is available at 
http://shoals.sam.usace.anny.mil/hawaii/pages/Oahu.htm). After investigating the origins 
of the dataset, however, and determining that these data were not collected by SHOALS 
in 2000, it is now believed that the data were collected by USGS circa 2002 to 2005. 
This information was provided by C. Fletcher (unpublished data, 2012). Prior to the 
completion of this particular project, there was no resolution concerning the actual details 
about the dataset. 

3. Software 

a. Environment for Visualizing Images + Interactive Data 
Language 4.8 (ENVI 4.8 + IDL) 

Environment for Visualizing Images (ENVI) 4.8 was the main software 
program used for data analysis and manipulation. ENVI is an image processing software 
system designed for multi- and hyperspectral data analysis and information extraction. 
This software is written in the Interactive Data Language (IDL), a programming language 
that provides integrated image processing and display capabilities [Research Systems, 
2004; Camacho, 2006]. 

ENVI was used to process the WV-2 imagery. All image mosaicking, 
spatial subsetting and chipping, radiance calibration, and registration was performed in 
ENVI. A number of IDL programs were written for each combination of two WV-2 
bands tested, and each incorporated a section to create masks for land, clouds, whitecaps, 
and glint, as well as apply a band ratio method for depth determination (see Section C of 
the Appendix). 

B. METHODS 

A number of pre-processing steps were perfonned before bathymetric derivation 
could occur. First, images were mosaicked and analyzed for overall quality and cloud 
cover. The next steps included map coordinate conversion, radiance calibration, 
land/cloud mask creation, sun glint removal, and the application of a ratio method for 
bathymetry derivation. 
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1. Bathymetry Derivation 

a. Assess Data Coverage, Angles, and Quality 

Before any data preparation steps occurred, images were assessed. It was 
important to ensure that imagery were collected over the correct site, and that datasets did 
not include any bad data. Scenes were ordered by mean satellite elevation angle (from 
most forward-looking, to most nadir, to most backward-looking). Table 3 and Section B 
of the Appendix list all 39 images acquired over Oahu. 
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Table 3. 


firstLineTime 

(Z) 

meanSatEI 

meanSatAz 

21:22:49 

24.9 

18.8 

21:22:58 

26.5 

19.2 

21:23:07 

28.1 

19.6 

21:23:16 

29.8 

20.1 

21:23:25 

31.7 

20.6 

21:23:34 

33.7 

21.2 

21:23:53 

38.2 

22.6 

21:24:02 

40.8 

23.5 

21:24:12 

43.7 

24.6 

21:24:21 

46.7 

26.0 

21:24:31 

50.1 

27.6 

21:24:41 

53.6 

29.7 

21:24:51 

57.6 

32.4 

21:25:01 

61.6 

36.1 

21:25:12 

66.0 

41.5 

21:25:22 

70.2 

49.2 

21:25:32 

74.0 

61.2 

21:25:42 

76.9 

80.1 

21:25:52 

77.8 

105.5 

21:26:02 

76.3 

129.8 

21:26:12 

73.1 

146.9 

21:26:23 

69.0 

157.9 

21:26:33 

64.8 

164.7 












b. 


Mosaic Data 


Data files were delivered in such a way that scenes needed to be 
mosaicked. This step was performed to be able to better view the entire coastline. 
Spatial subsetting could then occur without worrying about areas that may not have full 
coverage because the rest of the scene was in a different file. Figure 18 shows example 
images after mosaicking. 



Most Forward-Looking 
(24.9 degrees) 


Most Nadir 
(77.8 degrees) 


Most Backward-Looking 
(24.9 degrees) 


Figure 18. The most forward-looking, nadir, and backward-looking images 

after data coverage assessment and mosaicking 


c. Convert Map Projection 

Level IB (Basic) images were used, meaning that data were not projected 
to a plane using a map projection or datum and needed to be converted to Universal 
Transverse Mercator (UTM), World Geodetic System 1984 (WGS-84), Zone 4N 
(information is available at http://www.digitalglobe.com/downloads/ 
DigitalGlobe_Core_Imagery_Products_Guide.pdf). Pixels were also made square. 
Figure 19 shows how the map projection conversion affected the look of the imagery. 
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Figure 19. Map projection conversion from Geographic Latitude/Longitude (left) 

to UTM, WGS-84, Zone 4N (right) 


d. Radiance Calibration 

The WV-2 spectral radiance response is defined as the ratio of the number 
of photo-electrons measured by the system, to the spectral radiance [W-m -sr -pm ] at a 
certain wavelength present at the entrance to the telescope aperture. The spectral 
radiance response for each band is normalized by dividing by the maximum response 
value for that band to arrive at a relative spectral radiance response [Updike and Comp, 
2010 ]. 

Relative radiometric calibration and correction are necessary. This is 
because a uniform scene does not create a uniform image when it comes to raw digital 
numbers (DNs). This type of correction minimizes image artifacts, such as vertical 
streaks or bands due to differences in gain or offset, and is performed on raw data from 
all detectors in all bands during the early stages of WV-2 product generation. The 
products are linearly scaled to absolute spectral radiance [Updike and Comp, 2010]. 

In the case of large mosaics, radiometric balancing will help match the 
brightness of the other scenes used in the mosaic [Updike and Comp, 2010]. As glint 
removal is performed on images after the radiance calibration step, all images were left in 
radiance for this research. 
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e. Subset Data 

Images were subset to focus on Kailua Bay (Figure 20). This was done to 
reduce file size, emphasize a spot with interesting and variable bathymetry, as well as to 
better concentrate on a specific location that had less cloud cover. 



Figure 20. The full scene (left) was chipped to focus on Kailua Bay (right) 

/. Registration of Off-Nadir Images to the Most Nadir Image 

All non-nadir images were registered to the most nadir image. Twenty tie 
points were interactively selected for the same locations in all of the non-nadir images. 
The data were then warped using a first order polynomial. Maximum pixel error for the 
registrations was 1.521 pixels with an average root mean square (RMS) error of about 
0.895 pixels for the entire dataset. Images were then re-chipped so that every scene 
covered the same geographic coordinates. All chips were reduced to 995 samples and 
999 lines. An example is shown in Figure 21. 
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1010 


2010 


2100 


3100 



4100 


Figure 21. After registration, each image (top) was chipped to the 

995 x 999 pixel scene shown (bottom) 


g. Land, Glint, Cloud, and Whitecap Masks 

Application of spatial masks for land, glint, cloud, and whitecap areas 
allows analysis to be limited to only the areas and materials of interest. Masks were 
determined by comparing scatter plots using the Blue and NIR-1 bands. Water regions 
are those in the lower left portion of the scatter plot, as can be seen in Figure 22. 



Figure 22. The original image (left) has a corresponding scatter plot (middle); 

specific sections of the plot are highlighted to show how they correspond to the 
image (right)—red is mostly land, green is mostly whitecaps and sun glint, 

and blue is water 
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The scatter plot classified land, glint, clouds, and whitecaps as any pixels 
with values greater than the user-defined points chosen within the plot. Every other pixel 
would, therefore, be considered not land, glint, clouds, or whitecaps and would not be 
masked (Figure 23). Figure 24 illustrates how Figure 23’s user-defined region affects the 
imagery. 


2D Scatter Plot Used for Glint Removal 



6 8 10 12 14 16 

Blue Band 


Figure 23. 2D scatter plot of the Blue vs. NIR-1 bands; the user-defined outline 

region includes pixels containing water, and everything else can generally be 
classified as land, clouds, glint, whitecaps, etc. (e.g., not water) 
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Figure 24. Land, glint, clouds, and whitecaps have been masked in this image 

h. Glint Removal 

As discussed previously, if the ocean surface is perfectly flat, the 
reflection of the sun appears as a bright, small portion of the sea surface. If the ocean 
surface is not flat due to stronger winds, however, parts of the water surface further from 
the center of the sun glint pattern will be at the required orientation to reflect sun light to 
the viewer. Therefore, the sun’s reflection will cover a larger part of the ocean surface, 
and will be made up of many tiny highlights that each reflect from a particular point on 
the surface from the sun. The brightness of each pixel equals the total brightness of all 
highlights in that pixel and is, therefore, proportional to the fraction of the sea surface at 
the right slope [Kay et al., 2009]. 

If the sun and sensor are treated as a point source and detector for a given 
viewing geometry, there is only one facet slope and orientation that is consistent with 
specular reflection. The sun, in reality, has an angular diameter of 0.53 degrees, so there 
is a range of possible slopes that can reflect light from some part of the sun’s disc into the 
sensor [Kay et al., 2009]. 

Steps to remove glint were applied to the two WV-2 bands that were the 
focus of each version of code (Section C of the Appendix). The method used was based 
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on the method revised by Hedley et al. [2005] after Hochberg el al. [2003]. Work by 
Hedley et al. [2005] establishes the linear relationship between NIR and VIS bands using 
a linear regression based on a sample of the image pixels. Over areas with underlying 
spectral brightness, such as deep water, one or more regions with a range of sun glint are 
selected. For each VIS band, all selected pixels are included in a linear regression of NIR 
brightness (x-axis) against the VIS band brightness (y-axis). If the slope of this line for 
band i is b„ then all the pixels in the image can be deglinted in band i by applying the 
following equation: 

R i = Ri- b[(R\u{ - MinNiR), (2) 

which means: reduce the pixel value in band i (Ri) by the product of regression slope (/?,) 
and the difference between the pixel NIR value (Rnir) and the ambient NIR level 
(Min N i R ). R’i is the sun glint corrected pixel brightness in band i. Min N i R represents the 
NIR brightness of a pixel with zero sun glint, and can be estimated by the minimum NIR 
found in the regression sample or as the minimum NIR value found in the entire image 
[Hedley et al., 2005], Refer to Figure 25 for a graphical interpretation of the Hedley et al. 
[2005] method. 



Figure 25. Graphical interpretation of the Hedley et al. [2005] de-glinting method 

(From Hedley et al. [2005]) 
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The analysis incorporates the slope of the regression line. A simplified 
example of the modified equation used in the code (Section C of the Appendix) is as 
follows: 

Green = Green - Slope of Green Regression Line * (NIR-1 - MinNIR-1) (3) 

Yellow = Yellow - Slope of Yellow Regression Line * (NIR-2 - MinNIR-2) (4) 

There are, however, a number of differences between the Hedley el al. [2005] method 
and the one used for this research. This code utilized only two of the six possible WV-2 
VIS bands at a time. It also focused on the entire, global scene, rather than a small, local 
portion (only deep water, for example). It also incorporated masking, which was not 
used in the Hedley et al. [2005] research [Hedley el al., 2005] (refer to Section C of the 
Appendix). 


i. Band Ratio to Determine Relative Bathymetry 

A two-step process was used to derive bathymetry. Relative bathymetry 
was determined by performing a band ratio method, and then derived bathymetry values 
were obtained by regressing relative bathymetry values against verified depth data (to be 
further discussed in the upcoming sections). Camacho [2006] calculated relative 
bathymetry using the natural log transformed reflectance values of the deglinted 
reflectance image [Camacho, 2006]: 

ln( 1000* b 1 )/ln( 1000*b2) (5) 

The relative bathymetric values in this research were extracted using the following 
expression, with a slight adjustment to the constant value used by Camacho [2006]: 

ln( 100 * b 1 )/ln( 100 *b2) (6) 

where bl is the first band used, and b2 is the second. This is a modification of the 
equation used by Stumpfet al. [2003] [Stumpfet al., 2003]. 
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Comparison of Water Depths to LiDAR Bathymetry 


The bathymetry data of Kailua Bay were acquired in shapefile format. 
They were then converted to a raster in the proper map projection, geographically li nk ed 
to the image chip, and then clipped to the 995 by 999 pixels-size. These data were then 
used as the “true” data. 

k. Derived Bathymetry 

Derived depth values were calculated by regressing the relative depth 
values with the actual depth values collected by the “true” bathymetry. The chi-squared 
and correlation values were additional IDL code outputs that helped detennine the 
“goodness of fit.” A lower chi-squared value is indicative of a better fit. Correlation 
values are interpreted as percentages, so the highest value represents the best fit. 

2. Analyze Effects of Collection Geometry on Water Depth Derivation 

The effects of collection geometry on water depth derivation were analyzed by 
running all 39 images through the processing and analysis approach described above. 
Chi-squared and correlation values were recorded, and images of derived depth versus 
actual “true” depth plots were compared. Results will be discussed further in the 
following section. 


41 



THIS PAGE INTENTIONALLY LEFT BLANK 


42 



IV. RESULTS 


A total of 39 bathymetric maps were generated from the WV-2 data using the 
previously described band ratio method. The accuracies of these maps were then 
evaluated using “true” bathymetry data from the University of Hawaii at Manoa. The 
final WV-2 derived bathymetry was produced by regressing the relative bathymetry 
values against the actual “true” bathymetric measurement. The chi-squared and 
correlation values were then compared to analyze effects of multi-angle acquisition on 
depth derivation. 

A. BATHYMETRY FROM ENTIRE IMAGE 

A total of 39 bathymetric maps from a WV-2 multi-angle sequence over Kailua 
Bay, Oahu, Hawaii were generated. The following band combinations were analyzed: 

• Coastal & Blue, 

• Coastal & Green, 

• Blue & Green, and 

• Green & Yellow. 

Each band combination was analyzed using the band ratio method to obtain relative 
bathymetry. These values were then regressed against the “true” bathymetry data shown 
in Figure 26. 
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Figure 26. “True” bathymetry data; depth is represented by the Rainbow scale 

(shallow water is red and deep water is black) 

The results of the regression were the derived depth values. These were scaled to 
meters and plotted. Derived depth plots from Image IDs 1010 (most forward-looking), 
2010, 2100 (most nadir), 3100, and 4100 (most backward-looking) are shown in Figures 
27 through 31 for each band combination. All images have the same scale, a range of 
depths from 0 (red) to -18 (black) meters. 
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1010 



Coastal & Blue Coastal & Green Blue & Green Green & Yellow 

Figure 27. Derived depths for each band combination for Image 1010 

(most forward-looking) 


2010 



Coastal & Blue Coastal & Green Blue & Green Green & Yellow 

Figure 28. Derived depths for each band combination for Image 2010 


2100 



Coastal & Blue Coastal & Green Blue & Green Green & Yellow 

Figure 29. Derived depths for each band combination for Image 2100 

(most nadir) 
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3100 



Coastal & Blue Coastal & Green Blue & Green Green & Yellow 

Figure 30. Derived depths for each band combination for Image 3100 


4100 



Coastal & Blue Coastal & Green Blue & Green Green & Yellow 


Figure 31. Derived depths for each band combination for Image 4100 

(most backward-looking) 


It can be seen that none of these derived depth maps come close to a perfect 
correlation to the “true” bathymetry data. This most likely has to do with a number of 
factors, including high off-nadir acquisition angle and also cloud cover (seen as white, 
masked data in Figures 30 and 31). Due to the fact that wave patterns can be seen in the 
derived depth images, it is also speculated that a better glint removal method may be 
required in order to derive more accurate depths. 

Examination of these images shows that the Coastal & Green and Blue & Green 
WV-2 band combinations tend to perform better. Also, derived depth more closely 
matches the “true” depth as the sensor is acquiring images closer to nadir. 
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B. VARIATION IN ACQUISITION ANGLES 

In order to better understand the effect of acquisition angle on depth 
determination, plots of derived depth versus “true” depth for Image IDs 1010 (most 
forward-looking), 2010, 2100 (most nadir), 3100, and 4100 (most backward-looking) 
were created and compared. The chi-squared and correlation values of each were 
detennined and are displayed for Image ID 2100 (most nadir) in Figures 32 through 35 
(the rest can be located in Section D of the Appendix). 


2100: Coastal Sc Blue Derived vs. True Bathymetry 



-8 -7 -6 -5 -4 -3 -2 

Derived Bathymetry (m) 


Figure 32. Coastal & Blue: Derived depth vs. “true” depth for Image 2100 

(most nadir) 
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2100: Blue Sc Green Derived vs. True Bathymetry 



Derived Bathymetry (m) 


Figure 34. Blue & Green: Derived depth vs. “true” depth for Image 2100 

(most nadir) 
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2100: Green & Yellow Derived vs. True Bathymetry 



Derived Bathymetry (m) 


Figure 35. Green & Yellow: Derived depth vs. “true” depth for Image 2100 

(most nadir) 

Chi-squared and correlation values were calculated and plotted for each Image ID 
(view angle) for all four band combinations. These values represent the correlation of the 
depth derived from the band combination ratios versus the “true” depth (example values 
are shown in Section E of the Appendix). Results are depicted in Figures 36 through 39. 
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Coastal & Blue Chi-Squared 



Coastal & Blue Correlation 



Mean Satellite Elevation Angle 

Figure 36. Coastal & Blue: Chi-squared (top) and correlation (bottom) 

values (y-axis) plotted against the mean satellite elevation angle (x-axis) 

for all 39 WV-2 images 
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Coastal & Green Chi-Squared 



Coastal & Green Correlation 



Mean Satellite Elevation Angle 


Figure 37. Coastal & Green: Chi-squared (top) and correlation (bottom) values 
(y-axis) plotted against the mean satellite elevation angle (x-axis) for 

all 39 WV-2 images 
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Blue & Green Chi-Squared 



Blue & Green Correlation 



10% -I - 1 -i-i-.-i-,-.-,- 1 -i- 1 -i— 

24.9 29.8 382 46.7 57.6 702 77.8 69.0 56.6 46.5 38.7 32.5 27.4 

Mean Satellite Elevation Angle 


Figure 38. Blue & Green: Chi-squared (top) and correlation (bottom) values 
(y-axis) plotted against the mean satellite elevation angle (x-axis) for 

all 39 WV-2 images 
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4.40E+06 


Green & Yellow Chi-Squared 



24.9 29.8 38.2 46.7 57.6 70.2 77.8 69.0 56.6 46.5 38.7 32.5 27.4 


Mean Satellite Elevation Angle 
Green & Yellow Correlation 



Figure 39. Green & Yellow: Chi-squared (top) and correlation (bottom) values 
(y-axis) plotted against the mean satellite elevation angle (x-axis) for 

all 39 WV-2 images 

Although exact chi-squared and correlation values depend on the cutoff values 
used in the 2D scatter plots for each band combination (refer to Figure 23), the trend was 
the same after multiple runs. For all four band combinations, the images acquired closer 
to nadir had the lowest chi-squared values and the highest correlations. The optimal 
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values were found when running the Blue & Green code, with a correlation value around 
71%. The most erratic results came from the Green & Yellow code. This may be 
because the Yellow band did not penetrate as deeply into the water. 

Low chi-squared values are also seen from Image IDs 4050 to 4100. This is most 
likely due to the fact that clouds were entering the scene, being masked and, therefore, 
leaving very little water from which to collect spectral infonnation for depth derivation. 
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V. CONCLUSIONS 


A. SUMMARY OF OBSERVATIONS AND ANALYSIS 

Thirty-nine WV-2 images collected at multiple angles over a coastal, shallow 
water environment were used to analyze the effect of varying view angle on bathymetry 
derivation. Following initial data preparation steps, the entire dataset was analyzed to 
determine what effect, if any, acquisition angle has on nearshore depth retrieval. The 
data processing and analysis consisted of steps to remove sun glint, a band ratio method 
to determine relative depth, and a regression to find derived depth. 

Accuracy of depth retrieval did, in fact, demonstrate an association with image 
acquisition angle. Images acquired at more off-nadir view angles proved to have lower 
correlation values to the actual “true” bathymetry data. This was shown by the increase 
in correlation values for more nadir images within the dataset, followed by a decrease as 
images became more backward-looking. Likewise, chi-squared values decreased as the 
code approached the nadir-looking images, and increased again after that. 

Four band combinations were compared in an effort to detennine which WV-2 
bands might be best for determining depth in shallow, coastal environments. These 
included: Coastal & Blue, Coastal & Green, Blue & Green, and Green & Yellow. Of the 
four, the Blue & Green band combination perfonned best, with the highest correlation 
value of approximately 71% between “true” and derived depth. 

The results of this research showed that more accurate bathymetric depths will be 
derived from images that have been acquired closer to nadir viewing geometry. 
Additionally, when using imagery from the WV-2 sensor, a combination of the Blue and 
Green bands will penetrate into the water in such a way that a higher accuracy of depth 
information will be obtained. 

B. RECOMMENDATIONS AND FUTURE WORK 

There are a number of factors that may have negatively affected the data and 
analysis. While the scene over Kailua Bay was chosen because it was the least cloudy 
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beach out of the 39 images in the dataset, it still had cloud cover (especially between 
Image IDs 4050 and 4100). This led to fewer data points that could be analyzed for those 
scenes. It would also have been preferable to have used a more recent dataset for “true” 
bathymetry. This particular collection is thought to have been acquired between 2002 
and 2005, and in an environment as dynamic as a coastal shoreline, bathymetry could 
have changed significantly. In the future, it would be advantageous to find a recent 
bathymetric dataset, and then acquire WV-2 data over that geographic location. 

There are a number of possibilities for improvement and growth in future research 
related to this particular topic. It would be beneficial to apply the analysis approach to 
the entire dataset scene, rather than just the image chips. It would also be interesting to 
use a dataset with a wider range of satellite elevation angles (more off-nadir to more 
nadir and back). It might also be fascinating to apply this method to nearshore areas with 
black sand to see whether or not sand color affects the results of the algorithm. 
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APPENDIX: CHARTS AND GRAPHS 


A. EXAMPLE WORLDVIEW-2 METADATA FILE (FOR IMAGE 1010) 


version = "23.6"; 
generationTime = 
productOrderld = 
productCatalogld 
imageDescriptor = 
bandld = "Multi"; 


2011-08-01T01:59:55.000000 
_01_P001" ; 


"052564471010 


= "202001008362B200"; 
"BasiclB"; 


z; 


panSharpenAlgorithm = "None"; 

numRows = 3072; 

numColumns = 9216; 

productLevel = "LV1B"; 

productType = "Basic"; 

numberOfLooks = 1; 

radiometricLevel = "Corrected"; 

bitsPerPixel = 16; 

compressionType = "JPEG2000"; 

jpegProfileName = "nga_np je_pan_.nl"; 


BEGIN_GROUP 

= 

BAND_C 

ULLon 

= 

-157.86756401; 

ULLat 

= 

21.36508648; 

ULHAE 

= 

275.36; 

URLon 

= 

-157.49196245; 

URLat 

= 

21.69660890; 

URHAE 

= 

9.00; 

LRLon 

= 

-157.49051787; 

LRLat 

= 

21.56876101; 

LRHAE 

= 

11.00; 

LLLon 

= 

-157.87041130; 

LLLat 

= 

21.22582540; 

LLHAE 

= 

15.00; 

absCalFactor = 9.295654e-03; 


effectiveBandwidth = 4.730000e-02; 
TDILevel = 24; 


END_GROUP = BAND_C 


BEGIN_GROUP = 

BAND_B 

ULLon = 

-157.86756401; 

ULLat = 

21.36508648; 

ULHAE = 

275.36; 

URLon = 

-157.49196245; 

URLat = 

21.69660890; 

URHAE = 

9.00; 

LRLon = 

-157.49051787; 

LRLat = 

21.56876101; 

LRHAE = 

11.00; 

LLLon = 

-157.87041130; 

LLLat = 

21.22582540; 


63 




LLHAE = 15.00; 

absCalFactor = 1.783568e-02; 

effectiveBandwidth = 5.430000e-02; 

TDILevel = 10; 

END_GROUP = BAND_B 
BEGIN_GROUP = BAND_G 

ULLon = -157.86756401; 

ULLat = 21.36508648; 

ULHAE = 275.36; 

URLon = -157.49196245; 

URLat = 21.69660890; 

URHAE = 9.00; 

LRLon = -157.49051787; 

LRLat = 21.56876101; 

LRHAE = 11.00; 

LLLon = -157.87041130; 

LLLat = 21.22582540; 

LLHAE = 15.00; 

absCalFactor = 1.364197e-02; 
effectiveBandwidth = 6.300000e-02; 
TDILevel = 10; 

END_GROUP = BAND_G 
BEGIN_GROUP = BAND_Y 

ULLon = -157.86756401; 

ULLat = 21.36508648; 

ULHAE = 275.36; 

URLon = -157.49196245; 

URLat = 21.69660890; 

URHAE = 9.00; 

LRLon = -157.49051787; 

LRLat = 21.56876101; 

LRHAE = 11.00; 

LLLon = -157.87041130; 

LLLat = 21.22582540; 

LLHAE = 15.00; 

absCalFactor = 6.810718e-03; 
effectiveBandwidth = 3.740000e-02; 
TDILevel = 18; 

END_GROUP = BAND_Y 
BEGIN_GROUP = BAND_R 

ULLon = -157.86756401; 

ULLat = 21.36508648; 

ULHAE = 275.36; 

URLon = -157.49196245; 

URLat = 21.69660890; 

URHAE = 9.00; 

LRLon = -157.49051787; 

LRLat = 21.56876101; 

LRHAE = 11.00; 

LLLon = -157.87041130; 
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LLLat = 21.22582540; 

LLHAE = 15.00; 

absCalFactor = 1.851735e-02; 

effectiveBandwidth = 5.740000e-02; 

TDILevel = 6; 

END_GROUP = BAND_R 
BEGIN_GROUP = BAND_RE 

ULLon = -157.86756401; 

ULLat = 21.36508648; 

ULHAE = 275.36; 

URLon = -157.49196245; 

URLat = 21.69660890; 

URHAE = 9.00; 

LRLon = -157.49051787; 

LRLat = 21.56876101; 

LRHAE = 11.00; 

LLLon = -157.87041130; 

LLLat = 21.22582540; 

LLHAE = 15.00; 

absCalFactor = 6.063145e-03; 
effectiveBandwidth = 3.930000e-02; 
TDILevel = 18; 

END_GROUP = BAND_RE 
BEGIN_GROUP = BAND_N 

ULLon = -157.86756401; 

ULLat = 21.36508648; 

ULHAE = 275.36; 

URLon = -157.49196245; 

URLat = 21.69660890; 

URHAE = 9.00; 

LRLon = -157.49051787; 

LRLat = 21.56876101; 

LRHAE = 11.00; 

LLLon = -157.87041130; 

LLLat = 21.22582540; 

LLHAE = 15.00; 

absCalFactor = 2.050828e-02; 
effectiveBandwidth = 9.890000e-02; 
TDILevel = 6; 

END_GROUP = BAND_N 
BEGIN_GROUP = BAND_N2 

ULLon = -157.86756401; 

ULLat = 21.36508648; 

ULHAE = 275.36; 

URLon = -157.49196245; 

URLat = 21.69660890; 

URHAE = 9.00; 

LRLon = -157.49051787; 

LRLat = 21.56876101; 

LRHAE = 11.00; 
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LLLon = -157.87041130; 

LLLat = 21.22582540; 

LLHAE = 15.00; 

absCalFactor = 9.042234e-03; 

effectiveBandwidth = 9.960000e-02; 


TDILevel = 24; 

END_GROUP = BAND_N2 
outputFormat = "NITF21NCDRD"; 

BEGIN_GROUP = IMAGE_1 
s atId = "WV 02"; 
mode = "FullSwath"; 
scanDirection = "Forward"; 

Catld = "1030 0100 0DC522 00"; 

TLCTime = 2011-07-30T21:22:48.765050Z; 
numTLC = 2; 

TLCList = ( 

(0, 0.000000), 

(2055, 0.411000) ); 

firstLineTime = 2011-07-30T21:22:48.765050Z; 
avgLineRate = 5000.00; 
exposureDuration = 0.0002; 
minCollectedRowGSD = 7.002; 

maxCollectedRowGSD = 7.136; 

meanCollectedRowGSD = 7.069; 

minCollectedColGSD = 6.088; 

maxCollectedColGSD = 6.168; 

meanCollectedColGSD = 6.128; 


meanCollectedGSD = 


6.582; 

meanProductRowGSD 

= 

7.110; 

meanProductColGSD 

= 

6.162; 

meanProductGSD = 

6. 

619; 

rowUncertainty = 

294 

o 

LO 

colUncertainty = 

145 

• 97; 


minSunAz = 
maxSunAz = 
meanSunAz = 
minSunEl = 
maxSunEl = 
meanSunEl = 
minSatAz = 
maxSatAz = 


meanSatAz 


minSatEl = 
maxSatEl = 


meanSatEl 


96.0; 

96.4; 

96.2; 

72.3; 

72.3; 

72.3; 
18 . 8 ; 
18 . 8 ; 
18.8; 
2 4.9; 
25.0; 
24.9; 


52.7; 

52.9; 


minlnTrackViewAngle 
maxInTrackViewAngle 
meanlnTrackViewAngle = 52.8; 

minCrossTrackViewAngle = 14.3; 

maxCrossTrackViewAngle = 14.4; 
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meanCrossTrackViewAngle = 14.4; 

minOffNadirViewAngle = 54.1; 

= 54.1; 

= 54.1; 


cloudCover = 0.311; 
resamplingKernel = "CC"; 
positionKnowledgeSrc = "R"; 
attitudeKnowledgeSrc = "R"; 
revNumber = 9481; 

END_GROUP = IMAGE_1 
END; 


maxOffNadirViewAngle 
meanOffNadirViewAngle 
PNIIRS =1.3; 
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INFORMATION ABOUT THE 39 WORLDVIEW-2 IMAGES - IMAGE ID, ACQUISITION TIME, SENSOR 
ELEVATION, SENSOR AZIMUTH, OFF NADIR VIEW ANGLE, AZIMUTH, AND ZENITH 


Image 

ID 

firstLineTime 

(Z) 

meanSatEl 

meanSatAz 

meanOffNadirViewAngle 

Azimuth 

Zenith 

1010 

21:22:49 

24.9 

18.8 

54.1 

-161.2 

125.9 

1020 

21:22:58 

26.5 

19.2 

53.1 

-160.8 

126.9 

1030 

21:23:07 

28.1 

19.6 

52.0 

-160.4 

128.0 

1040 

21:23:16 

29.8 

20.1 

50.8 

-159.9 

129.2 

1050 

21:23:25 

31.7 

20.6 

49.5 

-159.4 

130.5 

1060 

21:23:34 

33.7 

21.2 

48.0 

-158.8 

132.0 

1080 

21:23:53 

38.2 

22.6 

44.6 

-157.4 

135.4 

1090 

21:24:02 

40.8 

23.5 

42.6 

-156.5 

137.4 

1100 

21:24:12 

43.7 

24.6 

40.3 

-155.4 

139.7 

2010 

21:24:21 

46.7 

26.0 

37.8 

-154.0 

142.2 

2020 

21:24:31 

50.1 

27.6 

35.1 

-152.4 

144.9 

2030 

21:24:41 

53.6 

29.7 

32.1 

-150.3 

147.9 

2040 

21:24:51 

57.6 

32.4 

28.7 

-147.6 

151.3 

2050 

21:25:01 

61.6 

36.1 

25.2 

-143.9 

154.8 

2060 

21:25:12 

66.0 

41.5 

21.4 

-138.5 

158.6 

2070 

21:25:22 

70.2 

49.2 

17.4 

-130.8 

162.3 

2080 

21:25:32 

74.0 

61.2 

14.3 

-118.8 

165.7 

2090 

21:25:42 

76.9 

80.1 

11.7 

-99.9 

168.3 

2100 

21:25:52 

77.8 

105.5 

10.9 

-74.5 

169.1 

3010 

21:26:02 

76.3 

129.8 

12.1 

-50.2 

167.9 

3020 

21:26:12 

73.1 

146.9 

15.0 

-33.1 

165.0 

3030 

21:26:23 

69.0 

157.9 

18.6 

-22.1 

161.4 




UrstLineTime 



meanSatEl 

meanSatAz 

21:26:33 

64.8 

164.7 

21:26:43 

60.6 

169.4 

21:26:53 

56.6 

172.8 

21:27:02 

52.9 

175.3 

21:27:12 

49.6 

177.2 

21:27:21 

46.5 

178.7 

21:27:30 

43.7 

179.9 

21:27:39 

41.1 

180.9 

21:27:47 

38.7 

181.7 

21:27:56 

36.5 

182.5 

21:28:04 

34.4 

183.1 

21:28:13 

32.5 

183.7 

21:28:21 

30.7 

184.2 

21:28:29 

29.0 

184.6 

21:28:38 

27.4 

185.0 

21:28:46 

25.9 

185.4 

21:28:54 

24.5 

185.7 







21 : 22:49 


meanSatEl 



21 : 28:29 








Angle 


Azimuth 
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21 : 22:49 



21 : 28:46 














c. 


IDL CODE FOR MOST HIGHLY CORRELATED BAND COMBINATION (BLUE & GREEN) 


;Final code for Krista Lee's September 2012 thesis 
;Code written by: Krista Lee, Angie Kim, Professor Olsen 

r 

;Blue vs. Green 

r 

;The purpose of this code is to take 39 image chips over Kailua Bay, Oahu, HI from 30 July 2011 and apply a 
deglinting process, 

;along with a band ratio method to determine water depth. The images, acquired at different angles, have 
alreaad been calibrated 

;for radiance, atmospherically corrected, and registered to the most nadir of the images. This code will 
then apply a simple 

;deglinting technique to remove the effects of sun glint, as well as use Stumpf et al.'s band ratio method 
for bathymetry derivation. 

;These depth values will be compared to actual lidar bathymetry data from the University of Hawaii. 

• kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

r 

;LOAD COLOR TABLE 

loadct, 13 ;Loads color table — Rainbow 

tvlct, r, g, b, /get 
r (255) = 255 
g (255) = 255 
b (255) = 255 
tvlct, r, g, b 

device, decomposed = 0, retain = 1 

•kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

r 

;DEFINE DIRECTORY 

dir = 'O:\Krista_Lee\Thesis\IDL\FINAL' 

cd, dir 

•kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

r 

;DEFINE FILES 

files = FILE_SEARCH ( '*_mosaic_wgs84_rad_resize*chip.dat' ) 

nfiles = n_elements ( files) 

print, string (nfiles ) + ' files found: ' 
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print, files 


fmask = '2100_mask_b.dat' ;Created specifically for this image chip in ENVI 

;Size of images 
samples = 995L 
lines = 999L 
bands = 8 

;FOR LOOP TO REPEAT PROCESS FOR EACH IMAGE 

for file_number = 0, nfiles - 1 do begin ;Loops through each elements in file array 

file = files(file_number) 

out_dir = 'Image_' + strmid(file, 0, 4) 

spawn, 'mkdir ' + out_dir ;This creates the folder for the WV-2 image 
cd, out_dir 

spawn, 'mkdir b_g' ;This creates the next folder down for the band combinations 

cd, 'b_g' , CURRENT = old_dir 

help, old_dir 

outdir = old_dir + '\b_g\' 

cd, dir ;Go back to the main directory where all the datafiles are, etc. 

;DEFINE OUTPUT FILES 

;Use input filename to create output filename 

filel9 = files(file_number) ;String containing name of the current file (filel9 = current file) 

shortfilename = strmid (filel9, 0, 4) 

band_ratio_vs_true_depth = shortfilename + '_band_ratio_vs_true_depth.tif' 
regression_file = shortfilename + '_regression.tif' 
final_depth_file = shortfilename + '_final_depth.tif' 

der_vs_true_scatter_plot_file = shortfilename + '_der_vs_true_scatter_plot.tif' 

• kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

r 

;OPEN LAND MASK 

mask = bytarr (samples, lines) 
openr, 1, fmask 
forrd, 1, mask 
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close, 1 


land_mask = mask 

mask = reform (mask, samples, lines) 
mask = rebin (mask, samples, lines, bands) 

datal = fltarr (samples, lines, bands) 
openr, 1, file19 
forrd, 1, datal 

close, 1 

;DISPLAY LAND-MASKED IMAGE 

window, 1, xsize = samples, ysize = lines, title = '(1) Water Only (land masked) - '+filel9 
datal = datal * mask 

tvscl, datal(*,*, [4, 2, 1] ) , true = 3, order = 1 

;MASK ADDITIONAL PROBLEM PIXELS (GLINT, CLOUDS, ETC.) 

;User defines water region 

window, 2, title = '(2) Scatter Plot - Blue vs. NIR-1 - ' + filel9 

ti = datal(*,*,l) ;Blue 
tj = datal(*,*, 6) ;NIR1 

order = sort(ti) & ti = ti(order) & tj = tj(order) ;Not strictly necessary, makes oplot easier 
nele = n_elements (order) & index = 10* lindgen (nele/10) ;Plots every 10th point 

plot, ti(index), tj(index), psym = 3, title = file, xtitle = 'Blue', ytitle = 'NIR-1' 

print, 'Click the upper-right-hand corner of the region that defines water — for masking purposes' 

cursor, x, y, /down, /data 
delta = 0.1 

xx = [ 0, x, x, 0, 0] ;For lines that create box when clicked 

yy = [ 0, 0, y , y, 0] 
plots, xx, yy 

glint_clouds_etc = where (datal(*,*, 1 ) gt x or datal (*,*,6) gt y) 
gee = glint_clouds_etc 

ngee = where (datal(*,*, 1 ) It x and datal (*,*,6) It y) ;ngee: NOT glint, clouds, etc. 
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help, gee 

datal= reform (datal, samples*lines, bands) 
datal(gee, *) =0 

datal = reform (datal, samples, lines, bands) 

gce_mask = bytarr (samples, lines) & gce_mask(*) = 1 & gce_mask(gee) = 0 
land_gce_mask = land_mask*gce_mask 

window, 3, xsize = samples, ysize = lines, title = '(3) Water Only (land and clouds masked) - ' + 
tvscl, datal (*,*, [4,2,1] ) , true = 3, order = 1 

• ■k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k'k-k-k-k 

r 

;DEGLINT USING HEDLEY APPROACH (UTILIZES SLOPE) 

;Shows scatter plot, best fit line, slop, and line along the length of the range 
window, 4, title = '(4) Scatter plot and fit for deglint (NIR-1 vs. Blue) - ' + filel9 
ti = datal(*,*,l) ;Blue 

tj = datal(*,*, 6) ;NIR1 

;Check for zeroes 

zeroes = where (ti le 0) 
nz = where (ti gt 0) 

ti = ti(nz) 
tj = tj(nz) 

order = sort(ti) & ti =ti(order) & tj = tj(order) ;Not strictly necessary, makes oplot easier 

nele = n_elements ( order) & index = 10 * lindgen ( nele/10) ;Plots every 10th point 

plot, tj(index), ti(index), psym = 3, /ynozero, xrange = [0, 2.0], xtitle = 'NIR-1', ytitle = 
degree= 1 

coeff_b = poly_fit(tj, ti, degree, chisq=chisql, yfit = yfitb) 
oplot, tj, yfitb 

slope_b = coeff_b(l) 

xyouts, 0.2, .2, string (slope_b) , /normal 

water_irl_min = min(tj) 
water_irl_max = max(tj) 
xmin = min(ti) 

plots, [water_irl_min, water_irl_max], [xmin, xmin] 

datal (*,*,1) = datal (*,*,1) - slope_b * (datal(*,*, 6) - water_irl_min) ;Blue and NIR-1 bands 

window, 5, title = '(5) Scatter plot and fit for deglint (NIR-1 vs. Green) - ' + filel9 
ti = datal (*,*, 2 ) ;Green 
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filel9 


'Blue' 


; NIR1 


tj = datal (*,*, 6) 

;Check for zeroes 

zeroes = where (ti le 0) 
nz = where (ti gt 0) 

ti = ti(nz) 
tj = tj(nz) 

order = sort(ti) & ti = ti(order) & tj = tj(order) ;Not strictly necessary, makes oplot easier 
nele = n_elements ( order) & index = 10 * lindgen (nele/10 ) ;Plots every 10th point 

plot, tj(index), ti(index), psym = 3, /ynozero, xrange = [0, 2.0], xtitle = 'NIR-1', ytitle = 'Green' 
degree = 1 

coeff_g = poly_fit ( tj, ti, degree, chisq=chisql, yfit = yfitg) 
oplot, tj, yfitg 

slope_g = coeff_g(l) 

xyouts, 0.2, .2, string (slope_g) , /normal 

water_irl_min = min(tj) 
water_irl_max = max(tj) 
xmin = min(ti) 

plots, [water_irl_min, water_irl_max], [xmin, xmin] 

datal (*,*,2) = datal (*,*,2) - slope_g * (datal(*,*, 6) - water_irl_min) ;Green and NIR-1 bands 

• • k-k'k-k-k-k'k-k-k-k'k-k-k-k'k-k-k-k'k-k'k-k'k-k-k-k'k-k-k-k'k-k'k-k'k-k'k-k'k-k-k-k'k-k'k'k'k-k'k-k-k'k-k'k-k'k 

r 

;DEGLINT IMAGERY BY SUBTRACTING IR WAVELENGTH REFLECTANCE 
imlb = datal 

• : k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

r 

;OBTAIN RELATIVE BATHYMETRY BY RATIO OF BANDS METHOD 

;Use expression: In(1000*bl)/In(1000*b2), where bl = Blue band, b2 = Green band 
indexl = where (imlb It 1) & imlb(indexl) = 1 
nfac = 100. ;Camacho used 1000. 

reldepth = alog (nfac*imlb(*,*, 1) ) /alog (nfac*imlb(*,*, 2) ) 

window, 6, xsize = samples, ysize = lines, title = '(6) Blue/Green Band Ratio - Bathymetry - ' + filel9 
rl = -1.1 ;min(truedepth) 

r2 = -1.0 ;max(truedepth) 

tv, bytscl (reldepth*land_gce_mask* (-1) , min=rl, max=r2), order = 1 

bar = bindgen(256) 

bar = reform(bar, 1, 256) 
bar = rebin (bar, 40, 256) 
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xO = 65 ;xO, yO: define location of color bar within the image 
yO = 30 

xl = xO + 40; + 1 ;xl, yl: define location of white box drawn around color bar 

yl = yO + 255; + 1 
tv, bar, xO, yO 

;Displays color bar on the image 

plot, bar, /nodata, /noerase, color = 0, xrange = [0, 1], yrange = [rl, r2] * (-1), $ 
xticks = 1, yticks = 4, yminor = 1, charsize = 1.0, xstyle = 1, ystyle = 1, /device, $ 
pos = [xO-1, yO-1, xl, yl], xtickname = [' '], yticklen = -0.05 

xyouts, xl+25, yO+128, 'Blue/Green Band Ratio', orient = 90, size = 1.0, color = 0, /device, align = 

0.5 


;Close unnecessary windows 

wdelete, 1 ; Closes Water Only (land masked) window 

; wdelete, 4 ;Closes NIR-1 vs. Green scatter plot 

; wdelete, 5 ; Closes NIR-2 vs. Yellow scatter plot 

• ■k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k'k 

r 

; OPEN GROUNDTRUTHED WATER DEPTH IMAGE 

truedepth_file = 'KailuaRaster_DepthMeters_resize_WGS84_2mpixels_995x999.dat' 

truedepth = dblarr (samples, lines) 
openr, 1, truedepth_file 
forrd, 1, truedepth 

close, 1 

truedepth = truedepth * (-1.0) ;Make depth values positive to match reldepth 

rl = -18 ;min(truedepth) 

r2 = 0 ;max(truedepth) 

bad = where (truedepth gt rl * (-1)) ;Gets rid of values deeper than -20 m (rl) - assume these are bad 

data 

nodata_mask = bytarr (samples, lines) & nodata_mask ( * ) = 1 & nodata_mask(bad) = 0 
land_gce_mask = land_gce_mask * nodata_mask 

;Place all masks on truedepth 
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(- 1 ) 


truedepth_masked = truedepth * land_gce_mask * 
truedepth = truedepth_masked 

goto, skipgt ;Skip ground truth image - - comment out this line if you want a window 
window, 7, xsize = samples, ysize = lines, title = ' ( 7 ) True Bathymetry - ' + truedepth_file 
tv, bytscl (truedepth, min = rl, max = r2), order = 1 

bar = bindgen(256) & bar = reform (bar, 1, 256) & bar = rebin (bar, 40, 256) 

xO = 65 ;x0, yO: define location of color bar within the image 

yO = 30 

xl = xO + 40; + 1 ;xl, yl: define location of white box drawn around color bar 
yl = yO + 255; + 1 
tv, bar, xO, yO 

;Displays color bar on the image 

plot, bar, /nodata, /noerase, color = 0, xrange = [0, 1], yrange = [rl, r2]*(-l), $ 
xticks = 1, yticks = 3, yminor = 1, charsize = 1.0, xstyle = 1, ystyle = 1, /device, $ 
pos = [xO-1, yO-1, xl, yl], xtickname = [' '], yticklen = -0.05 

xyouts, xl+25, yO+128, 'Depth (meters) ', orient = 90, size = 1.0, color = 0, /device, align = 0.5 
skipgt: print, 'skipping plot' 

• ■k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k'k 

r 

;COMPARE RELATIVE BATHYMETRY MEASURE TO GROUND TRUTH BY TAKING RATIO OF RELATIVE DEPTH TO TRUE DEPTH 
;Ignore areas where trudepth = 0 (truedepth data don't exist at those points) 

truedepth(bad) = 0 

good = where (truedepth ne 0) & help, good 
depth_ratio = fltarr (samples, lines) 
depth_ratio = reldepth/truedepth 
depth_ratio(bad) = 0 

rl = -0.5 ;min(depth_ratio(good)) 

r2 = 1.0 ;max(depth_ratio(good)) 

goto, skip_rat ;Not really useful 

window, 8, xsize = samples, ysize = lines, title = '(8) Blue/Green Band Ratio vs. True Bathymetry - ' + 
filel9 
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im = bytscl(-l * depth_ratio, min = rl, max = r2) 
index = where (land_gce_mask eq 0) 
im(index) = 255 
tv, im, order = 1 

bar = bindgen(256) 

bar = reform(bar, 1, 256) 
bar = rebin (bar, 40, 256) 

xO = 65 ;x0, yO: define location of color bar within the image 

yO = 30 

xl = xO + 40; + 1 ;xl, yl: define location of white box drawn around color bar 

yl = yO + 255; + 1 
tv, bar, xO, yO 

;Display color bar on the image 

plot, bar, /nodata, /noerase, color = 0, xrange = [0, 1], yrange = [rl, r2], $ 
xticks = 1, yticks = 2 , yminor = 1, charsize = 1 . 0 , xstyle = 1, ystyle = 1, /device, $ 
pos = [xO-1, yO-1, xl, yl], xtickname = [' , yticklen = -0.05 

xyouts, xl+25, yO+128, ' B/G Band Ratio vs. Depth (m) ', orient = 90, size = 1.0, color = 0, 

align = 0.5 

im = tvrdftrue = 1) 
im = reverse (im, 1 ) 
im = reverse (im, 3) 
im = reverse (im, 1 ) 

write_tiff, outdir + band_ratio_vs_true_depth, im 
skip_rat: print, 'Skipping ratio' 

• •kic'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

r 

;REGRESS RELDEPTH VS TRUEDEPTH TO FIND ABSOLUTE (DERIVED) DEPTH VALUE 

x = reldepth(good) 
y = truedepth(good) 

result = regress (x, y, SIGMA = sigmal, CONST = const, CHISQ = chisq, yfit = yfit, corr = corr) 
print, 'Equation of regression line (reldepth vs. truedepth): ' 

print, 'Regress Result: ' + string (result) + 'Regress Constant: ' + string (const) 
print, 'Chi-squared: ', chisq, ' Correlation: ', corr 


/device. 
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YFIT = yfit2 ) 


print, 'Valid regression?' 
degree = 1 

Result2 = POLY_FIT(x, y, degree, CHISQ = chisq2, SIGMA = polysig, 
print, 'Equation of regression line (reldepth vs. truedepth): ' 

print, 'Polyfit result: ' + string (result2) 
print, 'Chi-squared: ', chisq2 
print, 'Valid regression?' 
help, polysig, yfit, yfit2 

;Uncomment to see the scatter plot and best-fit regression line - works, just slowly, 
window, 9, xsize = 900, ysize = 900, title = '(9) Blue/Green Band Ratio vs. True Depth - ' 
x = reldepth(good) 
y = truedepth(good) 

good2 = where (y gt -12) 
xx = x(good2) 
yy = y(good2) 

nele = n_elements (xx) 

index = indgen (nele/10, /long) 

plo , xx(index), yy(index), psym = 3, xtitle = 'Blue/Green Band Ratio', ytitle='True Depth 
file, $ 

yrange = [-15, 0], xrange = [1.0, 1.2], ymargin = [5, 4] 
oplot, x(index), yfit(index) 


slope = result2 (1) 
intercept = result2 


xyouts, 

1.02, 

-0.5 

xyouts, 

1.02, 

-1.0 

xyouts, 

1.02, 

-1.5 

xyouts, 

1.02, 

-2.0 

im = tvrd ( 

) 


im = 255b 

- im 



tv, im 


( 0 ) 

'Chi-squared: ' + string (chisq) 

'Correlation: ' + string (corr, format = "(f8.3)") 
'Slope: ' + string (slope, format = "(f8.1)") 
'Intercept: ' + string (intercept, format = "(f8.1)") 


filel9 


(m)' , title 
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im = tvrd(true = 1) 
im = reverse (im, 1 ) 
im = reverse (im, 3) 
im = reverse (im, 1 ) 

write_tiff, outdir + regression_file, im 

• kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

r 

;USE REGRESSION RESULTS TO CORRECT RELDEPTH AND CREATE DERIVED DEPTH 
absdepth = reldepth * result (0) + const 

rl = -18.0 ;0.0 

r2 = 0.0 ;(—1.0)*min(absdepth);0.0 ;2.0 

window, 10, xsize = samples, ysize = lines, title = '(10) Derived Bathymetry - ' + filel9 ;ABSOLUTE IS 
THE INCORRECT TERM — DERIVED IS CORRECT 

im = bytscl (absdepth, min = rl, max = r2) 
index = where (land_gce_mask eq 0) 
im(index) = 255 

tv, im, order = 1 

bar = bindgen(256) 

bar = reform(bar, 1, 256) 
bar = rebin (bar, 40, 256) 

xO = 65 + 50 ;x0, yO: define location of color bar within the image 

yO = 30 

xl = xO + 40 ;+ 1 ;xl, yl: define location of white box drawn around color bar 

yl = yO + 255 ;6; + 1 

tv, bar, xO, yO 

;Display color bar on the image 

plot, bar, /nodata, /noerase, color = 0, xrange = [0, 1], yrange = [rl, r2], $ 

xticks = 1, yticks = 3, yminor = 1, charsize = 1.0, xstyle = 1, ystyle = 1, /device, $ 

pos = [xO-1, yO-1, xl, yl], xtickname = [' , yticklen = -0.05, $ 

ytitle = 'Blue/Green' 

xyouts, xl+25, yO+128, 'Derived Depth (meters) ', orient = 90, size = 1.0, color = 0, /device, align = 

0.5 

im = tvrd(true = 1) 
im = reverse (im, 1 ) 
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im = reverse (im, 3) 
im = reverse (im, 1 ) 

write_tiff, outdir + final_depth_file, im 

• ■k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k'k-k-k-k-k-k 

r 

x = absdepth(good) 
y = truedepth(good) 

result = regress (x, y, SIGMA = sigmal, CONST = const, CHISQ = chisq, yfit = yfit, corr = corr) 
print, 'Equation of regression line (absdepth vs. truedepth): ' 

print, 'Regress Result: ' + string (result) + 'Regress Constant: ' + string (const) 

print, 'Chi-squared: ', chisq, 'Correlation: ', corr 

print, 'Valid regression?' 

degree = 1 

Result2 = POLY_FIT ( x, y, degree, CHISQ = chisq2, SIGMA = polysig , YFIT = yfit2) 
print, 'Equation of regression line (absdepth vs. truedepth): ' 

print, 'Polyfit result: ' + string (result2) 

print, 'Chi-squared: ', chisq2 

print, 'Valid regression?' 
help, polysig, yfit, yfit2 

;Uncomment to see the scatter plot and best-fit regression line - works, just slowly. 

window, 11, xsize = 900, ysize = 900, title = '(11) Derived Bathymetry vs. True Bathymetry' 

x = absdepth(good) 

y = truedepth(good) 

nele = n_elements (x) 

index = indgen (nele/10, /long) 

plo , x(index), y(index), psym = 3, xtitle = 'Derived Bathymetry (m) ', ytitle='True Depth (m) ', 
file, $ 

yrange = [-15, 0], xrange =[-15, 0], ymargin = [5, 4] 
oplot, x(index), yfit(index) 

slope = result2 (1) 
intercept = result2 (0) 

xyouts, -14, -1.0, 'Chi-squared: ' + string (chisq) 

xyouts, -14, -2, 'Correlation: ' + string (corr, format = "(f8.3)") 


title 
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lm 


im = tvrd () 
im = 255b - 
tv, im 

im = tvrd (true = 1) 
im = reverse (im, 1 ) 
im = reverse (im, 3) 
im = reverse (im, 1 ) 

write_tiff, outdir + der_vs_true_scatter_plot_file, im 

stop 
endfor 

end 
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DERIVED DEPTH VS. TRUE DEPTH PLOTS AND LINES OF BEST FIT 


1010: Coostol & Blue Derived vs. True Bothymetry 


1 010: Coostol & Green Oerived vs. True Bothymetr y 


0.157 


Derived Bathymetry (m) 


Derived Bothymetry (m) 


1010: Blue 6c Green Derived vs. True Bothymetry 


1 010: Green Sc Yellow Oerived vs. True Bothymetry 


Derived Bothymetry (m) 


Derived Bothymetry (m) 











True Depth (m) True Depth (m) 



2010: Coostol & Blue Oerived vs. True Bothymetry 


Derived Bothymetry (m) 



2 010: Coo9tol & Green Derived vs. True Bothymetr y 


Derived Bothymetry (m) 



2010: Blue Sc Green Derived vs .True 8othymetry 


Derived Bothymetry (m) 



2010: Green & Yellow Oerived vs. True Bothymetry 


Corrtfotton: 0.2U 


•5.5 -5.0 -4.5 -4.0 

Derived Bothymetry (m) 
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3100: Coostol & Blue Oerived v3. True Bothymetry 


3 100: Coo3tol & Green Derived vs. True Bothymetr y 


CorrMboi* 


Derived Bothymetry (m) 


Derived Bothymetry (m) 


3100: Blue Sc Green Derived vs. True Bothymetry 


3100: Green A Yellow Oerived vs. True Bothymetry 


Derived Bothymetry (m) 


Derived Bothymetry (m) 



























4100: Blue & Green Derived vs. True Bathymetry 
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TABLES OF EXAMPLE CHI-SQUARED AND CORRELATION VALUES 
FOR EACH BAND COMBINATION 

1. Coastal & Blue (with best fits highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

4092798.6 

0.163 

1020 

26.5 

4081785.6 

0.169 

1030 

28.1 

4059720.3 

0.173 

1040 

29.8 

4043658.7 

0.19 

1050 

31.7 

3975405.5 

0.212 

1060 

33.7 

3878340.7 

0.254 

1080 

38.2 

3874253.7 

0.273 

1090 

40.8 

3740692.8 

0.319 

1100 

43.7 

3598962.8 

0.337 

2010 

46.7 

3520323.3 

0.372 

2020 

50.1 

3353757.3 

0.433 

2030 

53.6 

3377771.5 

0.425 

2040 

57.6 

3248184.9 

0.462 

2050 

61.6 

3216444 

0.472 

2060 

66.0 

3170698.1 

0.475 

2070 

70.2 

3187990.2 

0.48 

2080 

74.0 

3097967.8 

0.467 

2090 

76.9 

3107493.1 

0.483 

2100 

77.8 

3082345.9 

0.488 

3010 

76.3 

3105183.7 

0.489 

3020 

73.1 

3096221.6 

0.496 

3030 

69.0 

3119432.1 

0.493 

3040 

64.8 

3179182.1 

0.483 

3050 

60.6 

3215832.1 

0.475 

3060 

56.6 

3334628 

0.445 

3070 

52.9 

3375361 

0.434 

3080 

49.6 

3523816.2 

0.385 

3090 

46.5 

3585509.4 

0.359 

3100 

43.7 

3700396.6 

0.311 

4010 

41.1 

3768118.2 

0.28 

4020 

38.7 

3858040.5 

0.233 

4030 

36.5 

3893903 

0.212 

4040 

34.4 

3917829.2 

0.16 

4050 

32.5 

3875953.7 

0.139 
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4060 

30.7 

3757712.1 

0.127 

4070 

29.0 

3647051.5 

0.133 

4080 

27.4 

3520708.9 

0.113 

4090 

25.9 

3355284.3 

0.093 

4100 

24.5 

3020893 

0.038 


Coastal & Green (with best fits highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

3911460 

0.262 

1020 

26.5 

3888263.3 

0.271 

1030 

28.1 

3882215.3 

0.277 

1040 

29.8 

3847912.4 

0.288 

1050 

31.7 

3755520.8 

0.313 

1060 

33.7 

3628560.2 

0.354 

1080 

38.2 

3546724.9 

0.39 

1090 

40.8 

3317653.4 

0.452 

1100 

43.7 

3132645.7 

0.475 

2010 

46.7 

3004689.7 

0.516 

2020 

50.1 

2762942.9 

0.573 

2030 

53.6 

2698408.7 

0.591 

2040 

57.6 

2585662.4 

0.611 

2050 

61.6 

2493514.7 

0.632 

2060 

66.0 

2406520.2 

0.644 

2070 

70.2 

2380175.1 

0.65 

2080 

74.0 

2331197.5 

0.651 

2090 

76.9 

2262424.3 

0.66 

2100 

77.8 

2242532.2 

0.665 

3010 

76.3 

2294378.5 

0.665 

3020 

73.1 

2291861.7 

0.663 

3030 

69.0 

2360830.5 

0.653 

3040 

64.8 

2429096 

0.636 

3050 

60.6 

2569034.7 

0.615 

3060 

56.6 

2718923.5 

0.588 

3070 

52.9 

2771285.8 

0.577 

3080 

49.6 

3030260.5 

0.519 

3090 

46.5 

3148438 

0.487 

3100 

43.7 

3326883.7 

0.434 

4010 

41.1 

3239760.9 

0.455 
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4020 

38.7 

3554008.1 

0.357 

4030 

36.5 

3580085.4 

0.345 

4040 

34.4 

3690842.1 

0.284 

4050 

32.5 

3740913 

0.23 

4060 

30.7 

3641605.2 

0.213 

4070 

29.0 

3543715.8 

0.212 

4080 

27.4 

3412602.3 

0.203 

4090 

25.9 

3204708.3 

0.143 

4100 

24.5 

2973309 

0.115 


Blue & Green (with best fits highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

3837181.9 

0.294 

1020 

26.5 

3799491.4 

0.31 

1030 

28.1 

3826809.8 

0.295 

1040 

29.8 

3788817 

0.312 

1050 

31.7 

3674079.4 

0.34 

1060 

33.7 

3539442.7 

0.386 

1080 

38.2 

3449363.9 

0.419 

1090 

40.8 

3170962 

0.485 

1100 

43.7 

3013485.2 

0.512 

2010 

46.7 

2834575.7 

0.552 

2020 

50.1 

2629851 

0.601 

2030 

53.6 

2483369.3 

0.632 

2040 

57.6 

2372775.6 

0.652 

2050 

61.6 

2245770.1 

0.675 

2060 

66.0 

2157476.7 

0.689 

2070 

70.2 

2122878.3 

0.696 

2080 

74.0 

2099771.8 

0.698 

2090 

76.9 

2047377.9 

0.705 

2100 

77.8 

2019917.5 

0.707 

3010 

76.3 

2040766.2 

0.706 

3020 

73.1 

2064788.8 

0.703 

3030 

69.0 

2128616.8 

0.695 

3040 

64.8 

2209021 

0.682 

3050 

60.6 

2394974.9 

0.65 

3060 

56.6 

2537364.6 

0.626 

3070 

52.9 

2676606.9 

0.597 
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3080 

49.6 

2883049.7 

0.551 

3090 

46.5 

2994151.5 

0.521 

3100 

43.7 

3242079.3 

0.458 

4010 

41.1 

3074725.8 

0.499 

4020 

38.7 

3446563.2 

0.393 

4030 

36.5 

3455884.2 

0.391 

4040 

34.4 

3618164.1 

0.315 

4050 

32.5 

3676137.1 

0.268 

4060 

30.7 

3603561.2 

0.241 

4070 

29.0 

3515796.4 

0.23 

4080 

27.4 

3360508.4 

0.24 

4090 

25.9 

3203685.6 

0.232 

4100 

24.5 

2942699.1 

0.198 


Green & Yellow (with best fits highlighted) 


IMAGE 

ID 

meanSatEl 

Chi- 

Squared 

Correlation 

1010 

24.9 

4204121.8 

0.022 

1020 

26.5 

4199128.1 

0.056 

1030 

28.1 

4193524.9 

0.021 

1040 

29.8 

4195615.6 

0.008 

1050 

31.7 

4137792 

0.073 

1060 

33.7 

4088413.4 

0.121 

1080 

38.2 

4154217.8 

0.085 

1090 

40.8 

4059838.4 

0.161 

1100 

43.7 

4029843 

0.107 

2010 

46.7 

3892704 

0.212 

2020 

50.1 

3869278.8 

0.242 

2030 

53.6 

3803366.4 

0.278 

2040 

57.6 

3726302.1 

0.31 

2050 

61.6 

3732941.5 

0.309 

2060 

66.0 

3647327.9 

0.337 

2070 

70.2 

3558505.7 

0.374 

2080 

74.0 

3538108.6 

0.361 

2090 

76.9 

3526190.9 

0.36 

2100 

77.8 

3496716.8 

0.362 

3010 

76.3 

3537665.1 

0.371 

3020 

73.1 

3527937.5 

0.368 

3030 

69.0 

3792433.1 

0.267 
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3040 

64.8 

3710290.4 

0.322 

3050 

60.6 

4123412.3 

0.082 

3060 

56.6 

3981819.9 

0.21 

3070 

52.9 

3747042.4 

0.314 

3080 

49.6 

4140895 

0.012 

3090 

46.5 

4033670.6 

0.151 

3100 

43.7 

4077367.7 

0.07 

4010 

41.1 

3828814.7 

0.252 

4020 

38.7 

4043303 

0.095 

4030 

36.5 

4030102.2 

0.105 

4040 

34.4 

3891947.4 

0.18 

4050 

32.5 

3954476.1 

0.035 

4060 

30.7 

3821397.1 

0.045 

4070 

29.0 

3701158.5 

0.079 

4080 

27.4 

3565573 

0.068 

4090 

25.9 

3328975.3 

0.074 

4100 

24.5 

3061447.5 

0.064 
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